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Abstract: It is expected that thermalization may be described within gauge/gravity 
duality by considering time-dependent configurations on the gravity side of the corre- 
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we investigate a simple time-dependent toy model in which a mirror moves in the ra- 
dial direction of Anti-de Sitter space. For this configuration, we establish a procedure 
for calculating two-point functions of scalar fluctuations, based on a WKB approxima- 
tion. We test our method on two sample trajectories for the mirror, and find that the 
singularity structure of the two-point functions is in agreement with geometric optics. 
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1. Introduction 

Gauge/gravity duality has been extremely successful at describing strongly coupled 
systems by mapping them to weakly coupled gravity theories. In particular, a strongly 
coupled quantum field theory in d dimensions at finite temperature in thermal equi- 
librium can be described in terms of supergravity in a (d + l)-dimensional background 
AdS^+i of Anti-de Sitter black hole geometry [ffl, fj, [|. Recently, this technique has 
been proven to be particularly useful for describing properties of the strongly coupled 
quark-gluon plasma. 

For further progress in holographically describing the quark-gluon plasma as cre- 
ated in heavy-ion collisions, it is essential to consider thermalization, i.e. the relaxation 
from an initial non-equilibrium state to a final state in thermal equilibrium. Within 
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gauge/gravity duality, this is expected to be modeled by time-dependent geometries. 
Describing time-dependent geometries is a difficult task in general and often requires 
heavy use of numerics. Different avenues have been pursued so far towards this goal: 

One approach is to study the collision of gravitational shock waves in Anti-de Sitter 
space. Efforts along this line include @, ^ Q 0, [|. Different aspects of heavy ion 
collisions have been studied in shock wave collision models, such as early time dynamics 
0, ||, H, entropy production and critical conditions on thermalization || |?| |l0|, |TT 



More recently, the initial value problem for the non-linear Einstein equations has been 
solved numerically for planar gravitational shock waves ||. Despite the success of this 
approach, information other than the one-point function of the stress tensor remains 
very difficult to obtain due to the complexity of the metric resulting from the collision. 
For an initial attempt made in this direction, see [12| . 

An alternative approach is to consider a collapsing matter distribution in Anti-de 
Sitter space. It is natural to assume that in order to model far from equilibrium pro- 
cesses in quantum field theories and the relaxation into equilibrium, the supergravity 
picture should describe the dynamical process of black hole formation from some ini- 
tially regular spacetime. A mathematically clearly arranged setup is e. g. given by an 
infinitesimally thin but massive shell which collapses to form a stable black hole as a 
final state. Any process of black hole formation will lead from some initial problem 
to the propagation of bulk fields which have to satisfy the infalling wave boundary 
condition at the final equilibrium state. Before that, the boundary conditions for the 
supergravity fields will be time dependent and determined by the particular scenario. 

So far, the collapsing shell scenario has been considered in the quasi-static approxi- 
mation where the shell is considered to move adiabatically [TJ], , which simplifies 
the calculations significantly. In these approaches, a collapsing thin shell geometry is 
probed by a scalar field or a graviton field. The dual boundary two-point correlator 
is calculated as a function of the radial position of the shell and thus describes the 
thermalization process. In ||16|| , Hubeny, Liu and Ragamani proposed a bulk-cone sin- 
gularity conjecture, which states that the two-point correlator becomes singular when 
the two boundary points are connected by a bulk null geodesic. The application of 
this conjecture to a collapsing shell model establishes a connection between a distinct 
signature for the boundary observables and horizon formation in the bulk. 

A slightly different avenue in this context has been followed by |T7], ITS]. The 
authors of |T7| demonstrate analytically that a weak scalar perturbation collapses to 
form a black hole. The authors of |L8| investigate the gravitational collapse of energy 
injected into Anti-de Sitter space and the formation of an event horizon by considering 
the evolution of locally anisotropic metric perturbations initially located near the AdS 
boundary. At a late stage, the numerical results match with the analytical solutions 
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based on an asymptotic expansion [|19|, [2(J. An interesting link between the thin shell 



and perturbative approaches described above has recently been established in [21, 22 



p3fl . The evolutions of entanglement entropy in d = 2 and d = 3 have been studied in 
a Vaidya metric describing a collapsing shell in the formation of a black hole. In |23[ 
the equal-time two-point function, Wilson loop and entanglement entropy are explored 
in various dimensions - Signatures of time dependence for the chiral condensate have 
recently been investigated in [24| by embedding a probe D7 brane into the late time 
boost invariant background of Jl9| . 

In this paper we start a program of studying gravitational collapse of a thin mat- 
ter shell beyond the quasi-static approximation, thus generalizing | I5[| . For clarifying 
a number of technical details in the study of two-point functions for time-dependent 
gravity duals, in this paper we consider a simple toy model for a time-dependent ge- 
ometry consisting of a mirror moving in the radial direction of AdS space. We impose 
Dirichlet boundary conditions at the position of the mirror and calculate the two-point 
function of a scalar field in this geometry. In the special case of the mirror moving 
with constant velocity, scaling symmetry of AdS space is preserved. This allows us to 
solve for the complete set of eigenmodes of the scalar wave in terms of scaling variables, 
which are reminiscent of the late time scaling variable of Janik and Peschanski [|H| (see 
also the recent [p5|). We compute the two point function based on the eigenmode de- 
composition and find that the singularities of the two-point correlator are related to the 
physics of bouncing light ray between the moving mirror and the AdS boundary. Thus 
the singularity structure of the correlator is determined by a geometric optics picture. 
Our results generalize the static mirror case considered in |26[ to the time-dependent 
case. 

We explore the geometric optics limit in more detail with a WKB analysis, which 
enables us to reduce a general PDE to an ODE and leads to a prescription for calculating 
the two-point correlator for general radial trajectories of the mirror. The final formula 
for the correlator is expressed as a Mellin transform involving the ratio of incoming 
and outgoing waves for each component in the eigenmode decomposition. We test this 
correlator prescription using two sample trajectories of the mirror. In the first case of 
a mirror moving with constant velocity, we reproduce the geometric optics limit of the 
two-point correlator found previously using scaling variables. In the second case where 
a "mirror" moves along a spacelike geodesic, we find the singularities of the correlator 
are consistent with the bulk-cone singularities conjecture [IB|. 

The paper is organized as follows. In section 2 we solve the eigenmodes of the scalar 
wave equation in the presence of a mirror moving with constant velocity. In section 
3, we derive an explicit formula for the time-dependent, spatially integrated two-point 
correlator in terms of the eigenmodes. It takes the form of a Mellin transform. We test 
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this formula by considering its vacuum limit, where the mirror is absent. In section 4, 
based on a WKB analysis we establish a more general prescription for the two-point 
function which allows for an arbitrary mirror trajectory We test the prescription by 
reproducing the geometric optics limit of the correlator obtained previously. We also 
consider a "mirror" moving along a spacelike geodesic and find the resulting corre- 
lator consistent with the bulk cone singularities conjecture [OJ. Some details of the 
computation are moved to the appendices. 



2. Moving Mirror in AdS^+i 

In this section, we calculate the scalar two-point function for a special trajectory of 

the mirror which preserves scaling symmetry. We work with the AdS^+i metric in the 
Poincare patch 

R 2 

ds 2 = — (-dt 2 + dx 2 + dz 2 ) , (2.1) 
and start with the standard scalar wave equation 

' A(V^^W = 0. (2.2) 



The usual way of solving the equation of motion is to consider a specific Fourier compo- 
nent <j)(t, z) = e iut (j)((jj, z) to reduce the PDE to an ODE. In other words, the equation 
of motion is simplified by focusing on an eigenfunction of the operator d%. We know 
from quantum mechanics that this is possible, since the operator dt commutes with the 
Lagrangian operator 



c = 4=fA^rd v ) = ^ 2 (-d 2 + v 2 + a 2 ) - ^^d z . (2.3) 



An alternative explanation for the commutation of d t and C comes from the fact that 
time translation symmetry is an isometry of AdS^+i- In the presence of a moving 
mirror along the radial coordinate, the time translation symmetry will be broken by 
the Dirichlet boundary condition on the mirror. 

Defining t = ~, x = -|, £ = ~ , we can write L in terms of dimensionless coordinates 
i, x and z. In the following the bar will be suppressed and we should keep in mind 
that physical quantities are measured in units of the AdS radius R. Although the time 
translation symmetry is broken, scaling symmetry can be preserved for some special 
mirror trajectories. The scaling symmetry is generated by the operator 

C x = x^d^ = t d t + x V + zd z . (2.4) 
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We can verify explicitly that [£, C x ] = 0. For simplicity, we focus on solutions which 
depend on (t, z) only. Solving the eigenvalue equation, we obtain 

C X <P = A0 = v x/2 f(u) , (2.5) 

with v = tz, u = - and / being an arbitrary function of u. An obvious choice of the 
mirror trajectory that preserves the scaling symmetry is t/z — Uo, i.e. a mirror moving 
with constant velocity l/u . The general solution to ( |2.3| ) is given by 



/(«) = au^ (u - f (i^, l±i ; a + Izit) 



(2.6) 



The solution <p(u,v) = v^f{u) is analogous to the solution in momentum representa- 
tion, <p(t,z) = e lU)t 4>(u , z) . Here, <ft(u,v) describes the common eigenfunction of the 
Lagrangian £ and the scaling operator C x . We can view A as playing the role of the 
frequency u. Writing the solution in the more familiar coordinates (t, z), we have 

,/ x , rf- 1 / . \ i i — d f 1 d 1 + d , 3 — d I — t/z 
(f)(t,z) = Az—(t- z) x+ —F -—,^—;\ + 



2 ' 2 

d-i . .^.i-d^fl—d 1 + d „ 1 + d 1— t/z\ 

+ W—,—;-A+— J . (2.7) 

Near the boundary at z = 0, we find that the scalar field given by (|2.7| ) behaves as 

<f>{t, z)~t x + --- + z d t x ~ d + ■■■ (2.8) 

The two exponents correspond to non-normalizable and normalizable modes in AdS, 
respectively. In the flat limit z — > oo, i.e. in the deep interior of AdS, fl2.7| ) becomes 

<f>{t, z) ~ z^ (a(X) (t + z) x+1 ^r + b(\) (t - z) x+1 ^ , (2.9) 

with a(A) and b(X) two A-dependent constants. Obviously, (t + z) x+l ^ and (t — z) x+L ^ 
in fl2.9|) correspond to outgoing and incoming waves. Comparing (|2.9|) to the flat limit 
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of ingoing and outcoming contributions to the scalar wave in momentum representation, 

(j)(t, z) = zh iuJt H { j\uz) oc z ^ e Mt+z) ) 

2 

<f>{t, z) = zie^Hfiuz) oc z^e iuj(t ~ z) , (2.10) 

2 

we conclude that A should take the value A = ^- + iA with A an arbitrary real number. 
We now rewrite A as A = A' — such that A' is purely imaginary. In terms of A', 
(|2.7|) becomes 

.d-i, \ \t f 1 d 1 I d A . 1 — t/z 

<j>(t,z) = Az— W—,—; A' + 1;—^- 

d-i, , v fl — d l + d l — t/z\ 

+Bz—{t + z) x F[^—,——--\'+l- — -!—). {2.11) 



For simplicity, we drop the prime on A in the subsequent. 

The Dirichlet boundary condition at u = t/z = Uq fixes the ratio of A and B, 



B \u + lj V 2 '2' ' 2 J V 2 ' 2 ' '2 

(2.12) 



3. The Two-Point Correlator 
3.1 Derivation of the correlator 

In this section, we will use the solution of the scalar in the bulk to compute the cor- 
relation functions of the dual operator in the boundary field theory. We are interested 
in computing the correlation functions in coordinate space. We perform the "Fourier 
space" analysis for the transformation to A space in general d dimensions instead of 
setting d — 4. We will follow [f2| for the computation of the two point correlator 1 . The 
two point correlator of the operator dual to a massless scalar in AdS^+i is given by 

where S is the action of the scalar field. Using a regulator near the boundary z — e, 



1 Note a subtlety involved in this procedure is elaborated in [gTj, but this does not affect our result 
for the massless scalar. 
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the resultant action reads 2 

1 /" - - 



.V, =^ / dzdtd d - v x^g» v d^,,< 



I [ dtd d ~ l x ' 



*™(t) _ / ^dtd* 



(3.2) 



The second term vanishes by the bulk equation of motion. Furthermore, by the Dirichlet 
boundary condition <fi vanishes at the locus of the mirror z m (t) and we are left with 



1 f dtd^(t,z) d ^ z) 



yd-l 



(3.3) 



Note that our scalar wave in the bulk has no dependence on spatial x. The vertex that 
couples the source <f>(t,x), the boundary value of <f>(t,x,z), with the operator 0(t,x) 
simplifies to 

dtd d ~ 1 x<l){t,x)0(t,x)= \ dt(f)(t) / d d ~ l xO{t,x). (3.4) 



(3.5) 



The functional derivative of the action with respect to (j)(t) gives 
( / d d - l xO(t,x) [ tf'-VO^z')) -- 



lim S^MrVol. 



z d - x 8(j){t', z) 



with Vol = J d d ~ l x being the spatial volume, since we assume a spatially infinite 
mirror that does not break translational invariance in spatial dimensions. The causal 
nature of ( |3.5| ) will be specified in each explicit example later and discussed in section 4. 



Expanding ( 2.11]) near z = 0, the solution to ( |2.2|) takes the form 



(t, z) = K(d } A, A, B) t x+ ^ + ■■■ + L(d, A, A, B) t x -^z d + 



d+l 



(3.6) 



with 



K(d, A, A, B) 
L(d,X 7 A,B) 



1 + d 

i\— r{d) 



2 J r(i±*) 



l- d 
1\ — 



£H0 
r(¥) 



r(^ + A) r(^-A)' 
r(i + A) A | _r(i - A) B 



r(i^ + A) iw-a) 



(3.7) 



and the ■ • • denote terms of the form W(d, X)z d ' with integer dl ^ d. We do not write 



2 We will suppress the overall normalization from the SUGRA action in this paper. 
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down W(d, \)z d ' explicitly, since they do not contribute to the two-point correlation 
functions ( |3.5|) . Using the eigenfunction ( |2.6| ) as a basis set, we can express an arbitrary 
wave <p(t, z) as a superposition of 4>\(t, z), namely 

<f>(t, z) = j 4>x(t, z)g{\)d\ = J <f>° x (t)g(\)d\ + ■ ■ ■ + z d J <j>l(t)g(X)dX + ■■■ 

= <j) (t)--- + z d (j) d (t) + --- . (3.8) 

Here g(X) describes the weighting function for the component (p\(t,z) with eigenvalue 
A. In the following we will look for an explicit expression for the weighting function 
g{X), then using that result we will be able to write <p d {t) as a functional of <p°{t). This 
is the crucial step to determine the two point-correlator coming from inserting ( |3.8| ) in 
(|375| ), giving 

( j d d - l xO(t, x)0(t>, 0)) = -rf^|y , (3.9) 
where contact terms from • • • in (|3.8|) are excluded. 



At the moment we do not care about the causal nature of the correlator ( |3.9| ), but 
later we will discuss it in subsection \5.2\ The relevant explicit decomposition of 4>°(t) 
and <ft d (t) defined in ( |3~8D are obtained from (|3.6|) as 



/+ioo 
K(d,X,A,B) t x+ ^r g(\)d\, 
-ioo 

/+ioo 
L(d,\,A,B) t x ~^ g(X)d\. 
-ioo 



(3.10a) 
(3.10b) 



Defining 4>°{t) = we identify ( |3.10a| ) as the inverse Mellin transform. This 

observation allows us to invert ( 3.10a ) using the Mellin transform and obtain g(X) via 
the relation 



K{d,X,A,B) g(X) 



1 

2tH 



4>°(t)t x+ ^dt 



2ni 



<j>°(t')t'- x 2dt', (3.11) 



where in the intermediate step, a change of variable t — 1/t' is used. Plugging (p. 11 
to ( |3.10b|) , we obtain 



i r+ioo 

At) = ^ \ 

2vrz J_ ioo 



L(d,X,A,B) \ 
K(d,X,A,B) 



d+l 



t'-\tt')-^<i> (t') dt'dX. 



(3.12) 
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Using (|3.9| ), we end up with the correlator in the integral representation 



i»«0(*, .)0(f , 0)> = -- ^ 7 ,- A ; AB) x , (3.13) 



with A, A, 5) and A, A, £?) defined in Q and A, 5 satisfying ( gig ) 



Before we proceed to the evaluation of this expression, several comments on the 
correlator are in order: i) As is common in Minkowski signature, the correlator obtained 



from ( |3.13|) will depend on the specific wave we use in the bulk. In the mirror geometry, 
the causal structure of the correlator is in general complicated. In a certain limit, it 
should reduce to the retarded (B = 0) or advanced (A = 0) correlator. We will 
see later that the limiting correlator does agree with those obtained in momentum 
space representation; ii) We have derived the correlator for uq > 1, which is defined 
for t,f > 0. Actually most formulae are equally true for u < —1, corresponding 
to t, t' < 0. We will however, focus on the case uq > 1 for definiteness in what 
follows; iii) Our moving mirror does not introduce any dissipation to the background, 
the correlator should therefore be real; iv) The correlator should be finite. We note 
the possible divergent factor ^rrgr when d = 2, 4, 6 • • • . This potential pole should be 

cancelled out by the A- integral. 3 

3.2 Different limits of the correlator 



Before evaluating the integral ( 3.13|) , let us look at the limit B = and A = 0. These 



correspond to the incoming and outgoing waves, respectively. We should expect (|3.13|) 



to give retarded and advanced correlators from experience in the momentum space 
representation. We show in the following that this is also true in the A representation, 
which will serve as a nontrivial check of our prescription. For definiteness, we choose 

TV \ I l-\-d \ \ 

t,t' > 0. At B = 0, the A dependent part of the integrand simplifies to — — (4) . 

1 (AH — ^— ) \l / 

The poles in the complex A plane are at A = —n — with integer n > 0. The integral 
is only nonvanishing when the integration contour is closed counter-clockwise, i.e. when 
t > t' . Summing over residues, we obtain for the correlator defined in (|3.9|) 



S(f,d{t) -e(t t>) 2d r( ^ } «u) 
5<f>°(t>) [ \t-t') d+1 r(^)r(d)' ^- i4J 

For A = 0, the situation is quite similar. The A-dependent integrand simplifies to 
p|i-d ^ (f) A ; wr th the poles located at A = n + with integer n > 0. Physically we 



3 In those dimensions, the scalar wave contains a logarithmically divergent term near the boundary, 
which is encoded in the definition of the Hypergeometric function. 
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require d > 1, and as a result the integral is only nonvanishing when the integration 
contour is closed clockwise, i.e. when t < t' . Summing over the residues, we obtain 



5<p d (t) _ , . 2 d r(^) 



for the correlator. Let us now compare this result with the standard momentum and 
frequency representation, in which the incoming wave is given by 



which gives rise to 



z d ' 2 Hf(ujz) co > 
z d / 2 H\\-coz) co < 0, 

2 



(3.16) 



(3.17) 



Fourier transforming the above back to coordinate space, we obtain 



^ d it) QU ^ r (^+ 1 ) r(l-f)cos(7r(d + l)) 

500(f) 1 tj (t_t^+i r (l + f) 2 rf 7r ' 1 j 

After writing the cos-term as a product of two T-functions and using some relations 
between the T-functions, it can be shown that ( |3.18| ) and ( |3.14| ) are indeed identical, 



demonstrating the equivalence of using the A or the co representation. 

Keen readers may have noticed that our correlator vanishes for odd d. This is true 
in the domain of time we are interested in, i.e. for t ^ t' . In Appendix B, we present 
explicit examples of spatially integrated correlators for d = 3 and 4 starting with a 
general formula for unintegrated correlators in CFT 4 . These examples clearly display 
the subtle difference between the calculation in odd and even dimensions. 

A similar analysis shows that the advanced correlator obtained from the momentum 
representation also agrees with ( p. 15 ). This boosts our confidence in the A representa- 



tion of the correlator, and we will extract the time-dependent correlator in the moving 
mirror background using this representation. 

We focus on the UV part of the correlator following the work by Amado and Hoyos 
p6|. Specifically, these authors show that the UV part of the correlator, i.e. the part 



4 At t — t' , there is a non-analyticity resulting from the lightcone non-analyticity of the unintegrated 
correlators. 
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obtained by considering only frequencies with \u\ ^> 1, has an equivalent description in 
terms of geometric optics in AdS. Moreover, they found that the singularities of the cor- 
relator correspond to the time when the light ray bouncing between the AdS boundary 
and the mirror hits the boundary This is a special case of the bulk-cone singularities 
conjecture by Hubeny, Liu and Rangamani JjiJ. The latter, originally formulated in 
global AdS space, states that the singularities occur when the two boundary points are 
connected by a bulk null geodesic, i.e. by a light ray trajectory. In the Poincare patch, 
the light ray will not return to the boundary without being reflected at the mirror. 
It should be stressed that while the authors of |26] use a static mirror in AdS which 



introduces an explicit scale to the boundary CFT, our mirror moves in such a way that 
scale invariance of the boundary CFT is still preserved. Therefore, the UV limit of our 
case amounts to summing over all residues in the complex A plane, with |A| ^> 1. As 
we will see, the geometric picture is robust in our moving mirror geometry. 

The UV part of the correlator is evaluated in appendix A. As explained there in 
detail, the causal nature of the correlator is related to the chosen integration contour 
of A. In particular, in the UV limit, which amounts in particular to using the relations 
( A.3 ), the retarded correlator is given by 



(<<«-'* (t, x)0(f, o)> a = - «« - o D<0^ (£) " 1 1™^^) 



_i_ r?r((i-l) r b „\ 

xe 2 [a ' 



^(±1271-) ^ rUd + r+1) 

r=0 



d -d-i ( 7vb i 
Kt-t')-nT(U')—g - , (3.19) 



a d+1 ■ ■ ~ \ a 
with a = In and b = In |- . B n (x) are the Bernoulli polynomials and 

e i2 ™ = e^\ |c| < -. (3.20) 



is 



Note that the square bracket is only a function of c, which according to ( 3.20 
a periodic function with period 1. As c — > 0, the first term in the square bracket is 
singular while the second term is regular. The correlator diverges when b/a = n, with 
neN. Another interesting observation is that for odd dimensions (d = 1, 3, • • • ), the 
correlator vanishes identically due to the (^-dependent prefactor, which is a remnant of 
the behavior of the vacuum correlator. 

The periodic divergence of the correlator in the moving mirror geometry is consis- 
tent with the expectation from the geometric optics limit proposed in p6| : Suppose 
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Figure 1: The correlator contribution g (^) as denned in ( 3.19| ) versus ^ at d = 4. g (^) 
has a period of 4tt. In every odd interval {2nn, (2n + 1)tt), the correlator is positive, while in 
every even interval ((2n + 1)tt, (2n + 2)ir), the correlator flips sign. 




Figure 2: A schematic picture of a moving mirror in the AdS background, and the trajectory 
of a bouncing light ray. The correlator we obtain displays a singularity structure in agreement 
with this geometric optics picture. 



the mirror starts moving with velocity 1/uq at t = from the AdS boundary. If at tk 
the light ray reaches the boundary, it is easy to see that the next time the light ray 
hits the boundary will be at t^+i = ^zjtk, thus In tj j^- = ln^^j, which amounts to 
b = a. Indeed at b/a = n when b is an integer multiple of a, we observe singularities 
in the correlator. A plot of the correlator is included in Fig.[|. The structure of the 



correlator deserves some explanations: The function g f^>) has a period of 4tc, which is 

I i7r(d— 1) / f> \ _i_ iTt{d— 1) 

due to the oscillatory factor e a ta _c J = e a n . The period of the singularities 
in the correlator is 2n because there are two singularities in one period of g (— )■ From 
( 3.19p , one might expect the period to be 7r, however half of the singularities vanish as 
we sum over the plus and minus contributions in (|3.19|) . This is an artifact of choosing 
the retarded correlator. The Feynman correlator e.g. will pick either a plus or a minus 
sign, displaying the other half of the singularities explicitly. We have also included a 
view of a bouncing light ray in the presence of the moving mirror in Fig. This figure 
shows schematically that the time separation of the singularities increases with time as 
the mirror moves further and further away from the boundary. 

4. Moving Mirror in the Limit of Geometric Optics 

4.1 The WKB approximation and the limit of geometric optics 

As we have derived in the previous section, the two point correlator in the UV limit 
contains singularities when the times are related by a geometric optics path. In this 
section, we make this connection more precise. In particular, we develop a prescription 
for calculating the two-point function for general trajectories of the mirror. As a test 
of our prescription, we reproduce the result of the two-point function obtained in the 
previous section. 

The geometric optics limit can be described by a WKB solution for the scalar wave 
in the AdS background. Let us write the scalar wave as <p — Ae l8 ^ e , with A and 9 being 
the amplitude and the phase of the wave. The essence of the WKB approximation 
is that the phase 9 varies much faster than the amplitude A. We plug <fi = Ae %e l e to 
the equation of motion □</> = and perform a series expansion in e. The equation of 
motion for A and 9 are given by the leading order and next-to-leading order terms 5 

-e 2 + 9' 2 = 0, (4.1) 

-2A9 + 2A'9' + (-9 + 9") A - d —^A9' = , (4.2) 

z 

where the dot and prime denote derivatives with respect to t and z, respectively. The 
first equation can be solved for 9 = 9±(t±z). They will be used to eliminate the bracket 
in the second equation. The latter is solved by A = z^~ A±(t ± z). The positive and 
negative sign solutions have obvious identifications as outgoing and incoming waves. 
They form two linearly independent solutions to the EoM. Therefore, we split the 



5 We focus only on spatially homogeneous solutions and discard derivatives with respect to x. 



solution of Ocf) = into 



d-1 / i6,(t+z) »9_(t— *)\ d-1 

</>(t, z) = z— (A+(t + z)e — i — + A_(t - z)e — ; — J = z~ (F+(t + z) + F_(t - 2)) . 

(4.3) 

The WKB solution breaks down near the singularities of the original equation of motion 
□0 = 0. In the case of the AdS background, the only singularity is the AdS boundary 
z = 0. Near the AdS boundary, a general scalar wave can be written as a superposition 
of eigenmodes, 



<f>(t,z) 



g_(\)z—{t - z) x F 



l-d 1+d 



;A + 1 



l-t/z 



+g + (X)z—(t + z) x F 



\ / > / l-d 1 + A | 1; 1 - t/z 



dX 



(4.4) 



2 2' 7 2 

where g±(X) are the weighting functions for the incoming and outgoing components 



with eigenvalue A, respectively. In the UV limit |A| — > 00, ( |4.4| ) can be nicely matched 
with the WKB solution. Since lim 7 ^ OC) F(a, (3; 7; z) = 1, equation ( |4.4| ) simplifies to 



<f)(t, z)=z 



d-1 
2 



g _(X)(t - z ) x + g+(X)(t + zY 



dX 



(4.5) 



and we identify 



F+(t + z) = J g+(X)(t + z) x dX 
F-(t-z) = J g-(X)(t - z) x dX 



(4.6) 
(4.7) 



with F + and F_ defined in ( |4.3| ). 

We now aim for an expression for the correlator in a more general time-dependent 
situation, where the mirror follows a general, not necessarily scale invariant trajectory. 
Before proceeding to the derivation of such an expression, let us take a moment to look 
at the setting of the problem. Since we are interested in a time-dependent scalar wave 
in the bulk, which amounts to solving a 2D PDE (in the absence of x dependence), 
we need to specify both boundary conditions and initial conditions in order to have 
a unique solution. The boundary conditions are provided by the source at the AdS 
boundary and the Dirichlet boundary condition at the mirror. There is no further 
initial condition. 

The mirror with trajectory z = t/uo is special in the sense that it allows us to 
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study the "steady wave" (analogous to a planar wave in the AdS background), which 
does not require an initial condition. This particular trajectory does not couple the 
eigenmodes with different A and allows us to determine a ratio 9 g + ^ for every A. In 
terms of the boundary field theory, the missing initial condition is encoded in the state 
on which the operator 0(t, x) acts. 



In non-equilibrium field theory, the correlator should be studied using the in-in con- 
tour |28| , which is composed of one forward and one backward contour in the complex 
time plane. This was formulated holographically in real-time gauge/gravity duality by 
Skenderis and van Rees f29| . The state should be prepared by inserting sources in the 
Euclidean segments of the in-in contour. One of the essential points of |29| is that the 
bulk solution should be completely fixed by matching the Lorentzian segments to the 
Euclidean ones. The matching effectively provides the initial condition. Our mirror 
can be viewed as an effective source which creates the state. Despite the ambiguity 
related to the missing initial condition as described below ( |4.7| ), we will see that taking 
the UV limit as before actually allows us to obtain the wave solution without specifying 
an initial condition. In the UV limit |A| — > oo, the solution in the bulk only depends 
on t ± z, as defined in (Of), thus reducing the PDE to an ODE. 



To solve the wave equation, we start with the Dirichlet boundary condition on the 
trajectory of the mirror which implies that the wave vanishes there, 



<7_(A) (t - z m {t)) X + g+(\)(t + z m {t)Y 



d\ = 0. 



(4.8) 



Here and below, we suppress the integration bounds, with the understanding that the 
A integral runs always from — too to -Moo, while the t integral runs from to oo. This 
is to be combined with the asymptotic behavior near z = similar to ( |3.6p , where 
the unique weighting function g(X) in ( |3.8| ) is replaced by the more general weighting 
functions g±(X) as described in 



At) 



d-l 
1\ — 



r(^) 



(K + (A) </_(A) + K. (A) <7 + (A)) t x+ ^d\ , (4.9) 



with 



K + (X) 



r(i + A) 



A) 



r(i-A) 



E 



l+d 



A) 



(4.10) 



- 15 - 



For the sake of clarity, we introduce a shorthand notation 0(A) defined as 

0(A) = K+(X)g4X) + K4\)g + (\) , (4.11) 



hence equation Q4.9|) can be written as 

,0^ / M r ( rf ) 



u (t) 



r(^) 



0(A) t A+ ^rfA. 



(4.12) 



Next we will use equations (|4.8|) and the definition in (|4.11|) to write g±(A) in terms 
of 0(A). Later we will use the Mellin transform of (|4 .12 ) to express g±(X) as functions 
4>°(t). Then using these results, we can determine the coefficient function <p d (t) as an 
expression of 0°(t) and compute the two-point correlation functions. First using ( |4.11 
for rewriting (|4.8|), we obtain two relations 



J (K + (\)(t + z m ) x -K4\)(t-z m )") Yj^d\ = J 
J (K4\)(t-z m ) x -K + (\)(t + z m ) x ) jMjd\= J 



0(A) 
K_(X) 

m 

K + (X) 



(t + z m ) x d\, (4.13) 
(t-z m ) x d\. (4.14) 



Defining Vp(X,t) as the inverse of the time-dependent part on the LHS of ( |4.13| ) with 
respect to integration over t, i.e. 



K + (X')(t + z m (t)f - K4X')(t - z m (t)f = F(X', t), 

1 „,A - A', 



J dtV F (X,t)F(X',t) = 



(4.15) 



The insertion of the i's is to remind us that A and A' are purely imaginary. The solutions 
for fl4.13|) and ( |4.14|) can be formally written as 



9-W 
K_(X) 

9+W 



J dtdX' V F (\,t)^L( t + z m ) x ' , 
-J dtdX' p F (A,t)-|^l(t- 2m ) v 



(4.16) 



K + {X) 

respectively. The next step is to determine the near boundary coefficient function 4> d (t) 
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in terms of <fi°(t). Similarly to (4.9), <p d (t) is given by 



At) 



r(¥) 



(L + (A) </_(A) + L_ (A) g + (X)) t x ~^dX , (4.17) 



with 



T(i^ + A) ' 
After using the Mellin transform of (|4.12| ) 

<f>{X" - 



r(V-A) 



(4.18) 



dt'(f>°{t')t' 



2ni F(d) \2 

to replace 4>{X') by 0°(t') in ( |4.16|) and then inserting ( |4.16| ) into ([4.17|) , we obtain 



(4.19) 



At) 



2 d T(-d)T(±¥) 



2ni T(^' 



d\d\'dt" N(X,X',t",z m (t")) / dt' <j)°(t')- 

(tt') — 



4-1 J,0/j-/\ L 1 



X+i-X! 



(4.20) 



AAA. A'. :„, ) ee V F (X, t") [ L + (X)K_(X) {t '' + - L_(X)K + (X) (f " ] " 



K_(X') K ' ^ y ' K+(X' 

Using (|3.9| ), we end up with the correlator in the integral representation, 

-2 d dT(-d)T(^ 



d^xOfaxpit^O)) 



2ni T(^)T(d) 



t x t'- x ' 



dXdX'dt"N(X,X',t",z m (t"))- 

(tt>) — 

(4.21) 



This is our main result. The main difference between (|4.21|) and 



are the integra- 



tions over A' and t" in ( |4.21|) which encode the motion of the mirror for non-constant 
velocity. In the following we will see that ( 4.21 ) and ( 3.13 ) indeed give the same results 
if the mirror moves with constant velocity along the radial coordinate. For this pur- 
pose, we apply the result ( |4.21| ) to the case of a mirror moving with constant velocity 
t/z m = u . Then F(X,t) and its inverse Pp(A,i) defined in ( [4.15 ) are given by 



(4.22) 
(4.23) 



F(X,t)= \K + (X) 1 + - 

u 



*-(A)[l--| K 



V F {X,t) 



1 



,-x-i 



2mK + (X)(l + ^-K.(X)(l-i-y 
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Plugging the above expression to ( |4.21|) , 
The final integral of A is identical to ( |3.13| ) 
to 1, which is precisely the UV limit which 
of the correlator. As a result, we are bound 



we see the integrals of t" and A' are trivial, 
up to setting the hypergeometric function 
we used in appendix A for the evaluation 
to reproduce ( |3.19| ). 



4.2 The Correlator in the presence of a non-constantly moving "mirror" 

In this section, we will illustrate the procedure for obtaining the correlator in the 
presence of a mirror moving with non-constant velocity. We will see that the results 
are again consistent with expectations from geometric optics. The procedure derived 
in the previous subsection works for any trajectory of the mirror, but in practice this 
is hard to deal with for complicated trajectories, since it involves the inversion of the 
integral operator in order to obtain T>p(X,t) from F(X,t). In principle, this is possible 
to calculate by applying the Fredholm theory |31|| . However, this can be extremely 
complicated. We aim at finding some special trajectories which can lead to significant 
simplifications of the inversion procedure. We will see that it is possible at least in the 
UV limit. 

Suppose we propose the trajectory with t — z m = ^r—, or equivalently z m = 
\/t 2 — 1, defined for t > 1. This is a spacelike trajectory, thus cannot be associated 
with a physical mirror. We should understand it as providing a boundary condition in 
the bulk. Since this is a spacelike trajectory, any light emitted from the boundary into 
the bulk after t = 1 will not catch the "mirror" and no reflection is expected. There- 
fore, we should not see any singularities in the two point correlator according to the 
geometric optics picture. We will confirm this by explicitly computing the correlator 
in the UV limit. 

Note that F(X,t) = K+(X)(t + z m ) x - K_(\)(t- z m ) x as defined in flPUp and K+(X) = 
K_(— A). With the proposed trajectory, we have 

F(X, t) = K+(X)(t + z m ) x - K + (-X)(t + z m y x = -F(-X, t) . (4.24) 

In the limit A — > oo, K + (X) = „!^ A t+L ~ A 1 ^, where the argument is fixed by 

l (AH — 2~) 

|arg(A)| < 7r. Writing A = iA, we have the explicit expression for F(X,t) 

F(X,t) = sgn(A) |A|^2i sin ^|A|t + — ~ — J (4.25) 
= sgn(A)|A |— 2% cos (\A\t K — — - J , 
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with t = ln(t + z m ). To solve for / dtV F {\' ,t)F(X,t) = jS(^-), we look at the UV 
limit of the orthogonality relation of Bessel functions 



xJ u (£ > x)J l , (£'x) dx 



'0 

2 /*°° 77- 7 / 77- 77-7/ 77- 

- / cos(£r ) cos^'z -)dx = 6(£ - O . (4.26) 

Jo 



2 4' 



2 4' 



For this case we are able to identify the kernel defined in ( |4.16|) explicitly. The expres- 
sion above implies that V F (X,t) = sg "l^ cos (|A|f — ^ 1 j? \ +Zm . This formula 



4 / TTl t-\-z n 



[A|~2~~2i 

is only defined on the imaginary axis, but has a natural and useful extension in the 
complex plane, 



V F (\,t) = -^-(e xi ' 



1 ~\~ Z 



K + (X) J t + z, 



2m V K_(X) 
To further simplify (|4.20|) , we define 

F(X, t) = K + (X)(t + z m f + K4X)(t - z m f 

and express (t ± z m ) x in terms of F(X,t) and F(X,t), 

(t ± z m ) x = ±F(X,t) + F(X,t) 
K T (X) 2K + (X)K4X) ■ 

Inserting ( 4.29|) into ( 4.20 ), we obtain 



(4.27) 



(4.28) 



(4.29) 



i d (t) 2 d - 1 r(-d)r 



l+d\ 



5(j) (t') 2m T(^)T{d) 



dfdXdX' 



K^X')K + {X>) 



(4.30) 



P(X)V F (X, t")F{X\ t") + Q(X)V F (X, t")F(X', t"" 



with P(X) = L + {X)K„(X) + L_(A)7Y + (A) and Q(X) = L+(X)K„(X) - L_(A)iY+(A). 
Using the definition of T> F , the P(X) term can be simplified to 



2 d - l T(-d)T(^) 

2m f(0)f(dj J \K + {X) ' K-{X)J \t> 



L+W + ±M) (^dx 



(4.31) 



fl r(^)r\d) { 9{t ^ (t - t') d+1 + 6{t ' t] (f - t) d+1 . 
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This is just the average of the retarded and advanced correlators in the vacuum. We 
will not include this piece since we are only interested in a state (trajectory) dependent 
contribution to the correlator, which is given by the Q(X) term in ( |4.30|) . Let us rewrite 
(|4.30|) without the P(A)-term in the explicit form 



2 d ~ x T(-d)T 



¥3 



6(/>°(t') 



(2my r(±=^)r(d) 



dtdXdX 



N l + d W l + d 

, t X 2" f- A 2" 



-\t 



K+(\')K-{\') \K4X) K + (X) 



x (e A % + (A / ) + e- A '^_(A' 



— A't ; 



L + (X)K_(X) - K + (X)L„(X) 



(4.32) 



with K±, L± given by ( [4.10| ), ( 4.18|) . We perform this integral in appendix [C[ with 




2.5 



3.0 



Figure 3: The correlator (4.33) as a function of t for i' = 2 at d = 4. The green dotted and 
red dashed lines represent the first and second terms in the square bracket, respectively. The 
blue solid line is their sum. 



the result 
(/ d^x TO(t,x)O{t',0)} 



+ 



r(d) r(if^) L (tt'-i) 1 ^ 2r(if^)r(i + d) 

„ / ,3 + (i , t' - * 
x F l + d,— — ,2 + d; 



(# - 1) t) (tt'-ty+ d _ 



(4.33) 



for the correlator defined in (|3.9| ). This result is displayed in Fig. ||. Note that we 
have specified the correlator as a time-ordered one. This is because we focused on only 
one segment in the in- in contour. The causal nature follows from the principle of the 
real-time gauge/gravity duality [^9 . 
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We see that the d-dependence of the correlator is qualitatively similar to the con- 
stantly moving mirror case. It is again finite for even dimensions and zero for odd 
dimensions, a reminiscent of vacuum correlator. Turning to the singularities of the cor- 
relator, we see that the first term has singularities at t — t' — 1. The hypergeo metric 
function in the second term can be expressed in terms of elementary functions, but 
this is not necessary for our purposes. Singularities appear when the argument mEm 
becomes 0, 1 or oo, corresponding to t = t', t = 1 or t' = 1 respectively. The singular- 
ity at t = t' is just the usual lightcone one. The singularities at t = 1 and t' = 1 are 
closely related to the starting time of the mirror. No other singularity is found in the 
correlator, consistent with the geometric optics expectation. This example provides a 
further nontrivial check for our approach. 

We note that the relatively simple procedure leading to the two-point correlator 
is related to the SO(l, 1) symmetry preserved by the "mirror" trajectory t 2 — z 2 = 1. 
We expect that a method for solving the wave equation exists which is similar to 
that of Section 2, introducing variables which make the symmetry manifest. It will be 
interesting to see whether the two-point correlator following from such a method agrees 
with the result given by the general prescription, away from the singularities. 



5. Conclusion 

To summarize, we have solved the time- dependent problem of a spatially homogeneous 
scalar wave equation in AdS space in the presence of a radially moving mirror. First 
we considered a mirror with constant velocity which preserves scaling symmetry. We 
used the complete set of solutions to obtain the spatially integrated time-dependent 
two-point function for the CFT state defined by the moving mirror geometry. We 
found the result to be consistent with a geometric optics picture, in agreement with 



the bulk-cone singularities conjecture of [16 



Moreover, we have determined the precise relation between the geometric optics 
limit and the WKB approximation to the scalar wave solution. The established con- 
nection allows us to solve the scalar wave equation beyond scaling symmetry. As the 
main results of this work, we have established a formula for the two-point function 
corresponding to a moving mirror geometry with arbitrary trajectory, valid in the UV 
limit. We have tested this formula by reproducing the two-point function in the case 
of a mirror moving with constant velocity obtained before. We performed a further 
nontrivial check by considering a spacelike trajectory of the "mirror" (which is not 
reached by the light ray in this case). We found that the singularities of the result are 
again consistent with a geometric optics picture. 
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Our analysis provides essential technical tools for the study of two-point function in 
time-dependent geometries. These tools may now serve as a starting point for studying 
the behaviour of two-point functions during thermalization in strongly coupled field 
theories. We describe some immediate avenues to pursue in the subsequent. 

In particular, most immediately the prescription for the two-point function which 
we have obtained for a moving mirror in AdS space can be generalized to the case of a 
gravitational collapsing geometry. Specifically, this will involve solving the wave equa- 
tion in a thermal AdS background outside the collapsing shell. This simply amounts to 
a change of the lightcone coordinate from t ± z which we used in the present paper to 
£± J with f(z) the warp factor of the thermal AdS metric. Moreover, the Dirichlet 
boundary condition used here has to be replaced by a matching condition of a scalar 
outside and inside the shell. The key ingredient of matching the bulk WKB solution 
to the near-boundary eigenmode decomposition ( |4.4| ) will remain unchanged, since the 
metric of the collapsing shell is asymptotically AdS. We expect that this generalization 
will allow us to go beyond the quasi-static limit for the collapsing shell considered in 
]r5| and probe the early time dynamics in the thermalization of strongly coupled fields, 
an important ingredient in understanding the early thermalization of matter produced 



in heavy ion collisions |34 



A further generalization of our results to the boost-invariant background of is 
also possible: The idea is again to match the bulk WKB solution to a near-boundary so- 
lution, which is expressed as a superposition of the eigenmodes in pure AdS background. 
Due to the additional proper time r appearing in the boost invariant metric, both the 
WKB solution and the eigenmode set consistent with scaling symmetry will take a 
different form from the one discussed in this paper. This will give rise to a proper time 
dependent correlator with rapidity and transverse space coordinates integrated out, 



instead of a time dependent spatially integrated correlator. A direct application to |19 



with infalling boundary conditions of the wave imposed on the horizon will allow us to 



obtain the quasi-normal modes beyond the adiabatic approximation used in [[35 



Finally, let us compare our results for two-point functions of a massless scalar field 
with light-like separation to previous discussions of space-like separation, for instance 



those of |37], ^T|]. There it was pointed out that for large masses, or alternatively 
for large conformal dimension, the correlator is dominated by the contribution of the 
space-like geodesic joining the two points. Our result is that the pole structure of the 
correlator is determined by the geometric optics limit, i.e. by the trajectory of a light 
ray, is an analogous statement for the light-like case valid also for the massless case 
corresponding to small conformal dimension. We expect that our proposed study of 
correlators in thermalization geometries along the lines developed in this paper will 
also allow for a further study of the role of geodesies for correlators in these geometries. 
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A. Evaluation of (3.13) 

The core of the correlator is the following integral 



' L A fe + Jl^ B fe + r ( i 



2m J yr(\+ 1 -^) r(i^-A) J \v(x + ^) r(^-A) 

For defmiteness, from (J2.12) we take 




\\t-iA~ ^ 1 + d x n 1 — Uq, 



A=(u + — ; -A + 1; , (A.2) 



a„,l - d 1 + d x , , 1 - ito. 



B = -K - 1) A F(— , — ; A + 1; 

First of all, we note that the integrand (denoted as F(X,t,t')) has the property: 
F(X*,t,t') = F*(X,t,t'). Combined with the fact that the integration path is the 
imaginary axis, we see the integral is manifestly real, which is consistent with the re- 
ality condition of the correlator. We will use the residue theorem to evaluate ( |A.1[ ). 
The UV part of the correlator is given by the contribution from poles with large |A|. 
For this, we do an asymptotic expansion of the integrand. The following properties of 
T- functions and Hypergeometric functions are useful fl32"] : 

r(A + p) T{— X + p) sm7r(A — a) 7^°° 

A branch cut at the negative real axis is needed, with the argument of A fixed by 
|ar<7£| < 71. Using the above asymptotic behavior, we obtain the integrand ~ X d (jj) x 
as A — > 00. If t > t', we close the contour counter-clockwise and the integral receives 
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contribution from poles in the left half complex plane, while if t < t', we close the 
contour clockwise then the integral receives contribution from poles in the right half 
complex plane. 



The possible poles in the whole complex plane are poles of the Gamma function, Hy- 
pergeometric function and roots of 



r ( A+1 ). A + _ r (- A + i ; 



r(A + ^ 



■B = 0. 



(A.4) 



Note that F(a, /3;7, z) as a function of 7 has the same singularities as IX7) [32 , we 
can show that all the poles of the Gamma function and Hypergeometric function are 
removable. Thus we are only left with roots of ( A.4 ). Due to the non-algebraic nature 
of (|A.4|) , finding analytic expression of all the roots is not possible. However, we can 
deduce a general property of the roots: ( |A.4| ) can be equivalently written as 



r(A + i) 

r(A + ±±^) W-i 



u Q + 1\ 2 1 - d l + d 

r I — - — , — - — ; —A + 1; — - — / 



2 ' 2 ' 

r(-A + i) 



(A.5) 



r(-A + ^) W0-1 



u + l V j,A-d l + d 1 

f — - — , — - — ; A + 1; — 



Mo- 



or R(X) = R(-X) with 



R(X) 



r(A + i) fuo + i 



r(A + i±^ 



Uq 



n^^-A + i;^). (A.6) 



We note that R(X*) = R*(X). It is easy to show that if A is a root of flPQ , -A, A*, -A* 
are also roots. We plot the left hand side of (|A.4|) in the complex A plane, and find 
the zeros lie nearly equally spaced on the imaginary axis. Therefore, we conclude that 
the roots must be purely imaginary. Now let us determine the asymptotic form of the 
roots. In the limit A — > 00 (A — > 00), ( |A.4j ) has the following asymptotic expression 



A 2 



> + 1) A - (it - l) A e ±4 ^ 



(A.7) 



and the root is given by 



^ + 2k 
A = ±m— — — , 
ln M+l 

MQ — 1 



(A.8) 



with integer k > 0. Our approximate roots indeed are consistent with numerical plots 
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in the sense that they are symmetric with respect to the real axis and equally spaced. 
Furthermore we expect ( |A.8| ) to be more accurate when In ^-j- — > 0, i.e. uq — > oo. As 
the mirror moves more and more slowly, essentially all modes are effectively UV. 



The poles lie precisely along the integration contour of A. In order to obtain a well 
defined result, we have to deform the contour to circumvent the poles. The ambiguity 
associated with the detour corresponds to the different causal natures of the resulting 
correlator 6 . In practice, it is easy to calculate the advanced correlator, for which we 
shift the integration of A slightly to the left. Then all the UV poles lie to the right of 
the contour. The integration contour has to be closed counter-clockwise, which requires 
t' > t. In this way, we can avoid the branch cut on the negative real axis. The residue 
at each root is obtained with asymptotic expressions as 



r(A+i) 
r(A+i^) 



A 



r(-A+i) 



r(-A- 



l-d 
2 



B 



d_ 

dX 



r(A+i) A , r(-A+i) 



r(A+i±4) 



r(-A+i±^ 



B 



^iird 



t'J In ^±± 

u -l 



(A.9) 



We are happy to see the emergence of the factor (1 — e T47rc( ), which will precisely cancel 
the pole from — 
sum of residues 



the pole from ^1-)^ ■ Denote a = In an d b — In |-. The correlator is given by the 



d d - l xO{t,x)0{t'^)) A = -d6(t' - t) ^2% u — n d \l"J (ttT^ 

a r ( 



k=0 



-2k 



2 

(A.IO) 

evaluated at A = ±in 2 . The subscript "A" stands for the advanced correlator. We 
find that the sum over k can be expressed in terms of the Lerch transcendent function 
$(z, s, a) 



d d - l xO{t,x)0{t'^)) A = -dJ2 e (t' ~ t){tt')*e 



2a f I" 



x rM)r (¥) ^ d -i. 



TC-^)T(d) 



4 



(All) 



6 The ambiguity is familiar in the standard calculation of the vacuum correlator. It can be fixed by 
a prescription of the integration contour of the frequency 



As d — > integer, this reduces to 



( j d^xOit, x)0(t', 0)) A = -dJ2 Otf ~ *)(*0^ 



x e 



_ i-K{d— 1) / ;> 



(--c) 



d\ 



a TC-^)T(d) 



-c) d+1 (±i27r) 



r=0 



r!(rf + r + 1) 

(A.12) 



where the B n {x) are the Bernoulli polynomials |33|]. The constant c is defined as 
e «27rc _ e 1 -^- w j^] 1 | c | <- i 'Jq obtain the retarded correlator, we note a useful property 
of the integrand: F(X,t,t r ) = F(—X,t',t). This leads to the following relation between 
retarded and advanced correlators, 



( / d^xOit, x)0(t', 0)) A = ( / d^xOit', x)0(t, 0)) R . 



(A.13) 



B. The Spatially integrated correlator in d = 3 and d = 4 

We begin with the unintegrated correlator with conformal dimension d. The retarded 
correlator is given by pi 



We note the retarded correlator only has support at t = r — \x\. Integrated with 
j d d ~ l x for d = 3, we obtain 



6i 

7 V(-(t-«e) 2 + r 2 ) 2 (-(t + ie) 2 + r 2 ) 5 



rf 2 ^(t, £) = - ( ? ^ ^ ; ^ - , u ; ^ ; ^ |~ . (B.2) 



In the limit e — > 0, this vanishes identically. 

For d = 4, we have, apart from rational function of t and r, also the logarithmic 
function, which makes the ie prescription relevant 

The terms in • • • drop out as the limit e — > is taken. The logarithmic terms give rise 
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to a finite contribution 



J d 3 xG R {t,x) = -0(t)^ . (B.4) 



This agrees with —d ^ ^ when ( |3.14j ) is used in d = 4 dimension. 



C. Evaluation of (4.32) 



The core part of ( (4.32 ) is given by 

r(A + i)r(-A + 1) r(A + i)r(-A + 1) \ r(A' + ^)r(-A' + 





r(A' + i)r(-A' + i) 
_ vf r(-A' + i) . A ._ v 



r(-A + i) r(A + i) \ r(A' + ^) r(-A' + ±±^ 



(C.l) 



Note that t = ln(t + y/t 2 — 1) > 0. We introduce an upper cutoff T to regularize the t 
integral. With t integrated out, ( |C1[ ) takes the following form 



d\d\' 



r(A + i)r(-A + 1) r(A + i)r(-A + 1) \ r(A' + ^)r(-A' + 



r(A + if2)r(ip _ a) r(A + i±^)r(if^ - A) y r(A' + i)r(-A' + 1) 



r(-A + i±^)r(A' + 1) e ( A + A ') T - i r(-A + ^)r(-A' + 1) e ( A - A ') T - 1 
r(-A + i)r(A' + ^) a + A' + r(-A + i)r(-A' + ^) a -A' 
_ r(A + ^)r(A / + i) e (- A+A, ) r - i _ r(A + ^)r(-A / + i) e (- A - A/ ) r -i 
r(A + i)r(A' + i±^) -a + A' r(A + i)r(-A' + W -A -A' 



(C.2) 



We now use the residue theorem to evaluate the integrals of A and A'. Note that the 
appearance of ^^llv" 1 ^ oes no ^ introduce any new poles as they have finite limits 
when ±A ± A' — > 0. We begin with the integral of A'. It is helpful to keep in mind 
that t, t' > 1 since the "mirror" does not leave the boundary until t = 1. Completing 
the A' integral, we find only the first and third terms in the bracket contribute and the 
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dependence on the cutoff T drops out naturally. We obtain the result 



2ni J dX 

-1) 



r(A + i)r(-A + i) 



r(A + i)r(-A + i) 



r(A + ±=^r(-A + m r(A + ^)r(-A + ^) 



X 



E 

n=0 



n! r(-n 



l+d 



+ 1) 



-t' 



l + d 



r(-A + ^ 



A + n+i±^ r(-A + l) 



+ 



(C.3) 

r(A + ^ 



-A + n + i±^ r(A + i) J 



We again perform the contour integral. This time it is much simpler: As t > 1, only 
the poles in the left half plane contribute. Considering the first term in the sum, we 
find that poles at A + 1 = — m for m e {0, 1, 2 • • • } and A + n + = are relevant. 
The contribution from the first set of poles is proportional to 



E 

m=0 



-iy 



ml 



T(m + 1 + 



r(-m-i + i^) r(m + i + i^)r(-m-i + i^; 



-m— 1 



m + 1 - n - 



sin 



M3±d)\ r ( 3+d 



mm 



71 



r(^) 



3 + d l-d 3-d 1 

— ^ — > ~ n ^ — ^ — j _n ~! — ^ 



(C.4) 



We note that the d-dependent prefactor vanishes identically by the properties of the 
Gamma function. The pole at A + n + = gives the contribution 



E 

n=0 



-l) n / 1 r(l + d + n) 

nF \T(-n - d) ~ r(l + n)T(-n) 



(tt'y 



{tt') — 



T(-d) (tt 1 - l) 1+a 



where 



(C.5) 

cancels the overall divergent factor r(— d) to yield a finite numerical coef- 
ficient. Next we consider the second term in the sum of (|C.3| ), the relevant poles are at 
A + = — m, with the contributions 



T(-d) 



E 

n,m= 
oo 

E 

n,m=0 



:-i) 



m+n 



r(m + ±^ + i) 



o m!n! r(-n - + l)T(-m - d)T(m + 1 + d) m + n + l + d 
-l) m+n sin vr(n + ^)T(n + ±±^) sin vr(m + d)T(m + + 1) 



mini 



71 



71 



t'- n t~ 
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(tt')- — 



n,m=0 



l + d 



(2 + (i) m+n m!n! 
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where (a) m = is the Pochhammer symbol. The double sum can be expressed in 

terms of Appell function, which again can be converted to a Hyp ergeo metric function 



sin7r^± smndTi^rm^) 



7T^ 



1 + d 



Fi l + d. 



3 + d l + d 



sin7r^isin7rrfr(i±^)r(^^ 



7T^ 



l + d 



i-l 

f 



2 

-i-d 



l + d 



(C.7) 



F(l + ^,2 + ,;^ 



(tty 



l+d 



Writting sin-functions as product of two T-functions, collecting all nonvanishing terms 
and inserting the overall factor, we end up with the following result 



S(j) d (t) 



ifl(t>) T(^)T(d) 
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2T(^)T(l + d) 
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